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ABSTRACT 

We report the discovery of the first low-mass planet to emerge from the NASA-UC Eta-Earth 
Program, a super-Earth orbiting the KO dwarf HD 7924. Keplerian modeling of precise Doppler radial 
velocities reveals a planet with minimum mass Mp sin i = 9.26 M® in a P — 5.398 d orbit. Based on 
Keck-HIRES measurements from 2001 to 2008, the planet is robustly detected with an estimated false 
alarm probability of less than 0.001. Photometric observations using the Automated Photometric 
Telescopes at Fairborn Observatory show that HD 7924 is photometrically constant over the radial 
velocity period to 0.19mmag, supporting the existence of the planetary companion. No transits were 
detected down to a photometric limit of ~0.5mmag, eliminating transiting planets with a variety 
of compositions. HD 7924b is one of only eight planets known with Mpsini<10M® and as such 
is a member of an emerging family of low-mass planets that together constrain theories of planet 
formation. 

Subject headings: planetary systems — stars: individual (HD 7924) — techniques: radial velocity 



1. INTRODUCTION 

The NASA-UC Eta-Earth Survey by the California 
Planet Search (CPS) group is a systematic search for low- 
mass planets (~3-30M®) orbiting the nearest 230 GKM 
stars suitable for high-precision Doppler observations at 
the Keck Observatory. In order to place statistically sig- 
nificant constraints on the fraction of stars with Earth- 
mass planets in the habitable zone, ?7®, and theories of 
planet formation, each star is observed a minimum of 20 
times at ~lms _1 precision, including at least one set of 
high-cadence observations on consecutive nights. Planet 
candidates are followed up with additional high-cadence 
observations to confirm the candidate signal and to ac- 
curately measure the planet's orbital parameters. While 
Doppler surveys cannot detect 1 M® planets in the hab- 
itable zone, we will estimate 77® by modest extrapola- 
tion from the distribution of detected super-Earths and 
Neptune-mass planets in shorter-period orbits. This pa- 
per is the first in a series to describe the planets and plan- 
etary systems emerging from the Eta-Earth Program. 

The statistics of planet occurance from the Eta- 
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Earth Survey will offer important const raints on 
comp e ting theories of p l anet fo rmation (llda fc Linl 
1 20041: iKenvon fc Bromlevl 120061: lAlibert et alj 120051: 
iMordasini et all l2008t llda fc Lini |200j). These models 
differ in assumptions about the growth rate of dust into 
planetesimals, the viscosity of the disk, the location and 
effects of the snow line in the disk, the efficacy of inward 
migration, the accretion of gas and water, and the rele- 
vance of planet-planet interactions. They are consistent 
with detections and mea s urements of jovia n gas giants 
(e.g. iMarcv et alj l2005at lUdrv et alj l2003h . However, 
these theories predict that planets of mass 1-30 M® are 
rare within 1 AU, forming a "planet desert" of super- 
Earths (Mpisini < 10 M®) and Neptune-mass plan- 
ets. All models predict that Type I migration quickly 
(r ~ 10 5 yr) causes the rocky planets to spiral inward, 
destined to be lost in the star 10 . Meanwhile, the more 
massive rocky cores (M > 1 M®) accrete gas quickly, be- 
coming ice- or gas-giants. The resulting distribution of 
planetary systems is depleted of planets in the mass range 
1-30 M® within 1 AU. 

If the prediction of a low-mass desert is contradicted 
by a statistically well-defined sample of stars (such as 
the Eta-Earth Survey), then planet formation theory 
must be significantly modified with new physics. The 
planet desert has indeed been challenged bv lMavor et alJ 
(12009ft who estimate that 30% ±10% of GK dwarfs 
have rocky or Neptune-mass planets inward of 50 d or- 
bits. This claim is based in part on detections by 
the Swiss group of three N eptune-size planets orbiting 
HD 69830 ijLovis et al.ll2006l) . two of three planets orbit- 
ing GJ581 in or near the Habitable Zone ()Bonfils et al.l 
l2005HUdrv et al.ll2007lh and a triple super-Earth system 
around HD 40307 (jMavor et alJl2009T >. 

It is against this backdrop of recent discoveries and 



10 While some of these models (e.g. Ida & Lin 2008) show a pile- 
up of low-mass planets in short-period orbits (analogous to hot 
Jupiters), this concentration may be artificial since the physical 
mechanism for stopping the migration remains uncertain. 
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TABLE 1 
Stellar Properties of HD 7924 



Parameter 



Value 



Spectral type 

M v 

B-V 

V 

Distance (pc) 
[Fe/H] 
T e g (K) 
^sini (kms ) 
log g 
U (L e ) 

m* (M e ) 
R* (Re) 

logij' K 

P rot (days) 
o-phot (mag) 



KOV 

6.056 

0.826 

7.185 
16.8±0.1 
-0.15 ±0.03 
5177 ±44 
1.35±0.5 

0.36Btffi 

n 000+0.022 
u.ooz_ Q 036 

0.78 ±0.02 

0.20 

-4.89 

~38 

< 0.0017 



competing claims of the abundance of rocky planets that 
we announce the detection of a super-Earth orbiting the 
K dwarf HD 7924. §[2] of this paper describes the host 
star, HD 7924. The spectroscopic observations and their 
Doppler reduction are described in §[3] Our detection 
of HD 7924b, including Keplerian fitting and false alarm 
probability estimation, is described in §|4] In §[5j we con- 
sider and ultimately reject a second planetary companion 
based on the current data. Our photometric observations 
and their constraints on planetary transits are described 
in §[6] We conclude in §[7] with a discussion of HD 7924b 
and the emerging family of super-Earth planets. 

2. PROPERTIES OF HD 7924 

HD7924 (HIP 6379, GJ56.5) is a KO dwarf whose 
properties are summarized in Table [T] It is 
nearby, d = 1 6 . 8 pc , an d relatively bright, V — 7.185 
(jPerrvman et al.l 119971 : Ivan Leeuwenl 120071 ). With 
My =6.056 and B - V = 0.826, HD7924 star lies 
0.27 mag b elow the Hipp arcos average main sequence as 
defined bv I Wrightl (|2005h . 

Using the SME (Sp ectroscopy Made Easy) LTE 
spectral synthesis code, IValenti fc Fischer] |2005) ana- 
lyzed a high-resolution, iodine-free Keck-HIRES spec- 
trum of HD 7924 and found the effective temperature 
T off = 5177±30K, surface gravity log£ = 4.58 ±0.08, 
projected rotational velocity v sinz = 1.35 ± 0.5 kms -1 , 
and iron abundance ratio [Fe/H] = —0.15 ± 0.03 (i.e. 
slightly metal poor, consistent with its location be- 
low the average main sequence). Using these quanti- 
ties and catalog values for U-band photometry, paral- 
lax, etc., they interpolated the Yonsei-Yale isochrones 
to obtain the luminosity L+ = 0.368iooi,8 L®, radius 
R* = 0.754 ± 0.013 R Q , and mass M* = 0.80tj^ M . 
Isochrones provide no useful constraint on the age be- 
cause the star has not begun evolving off the main se- 
quence; 

iTakeda et al.l (|2007t ) used th e spectroscopically deter- 
mined stellar parameters from IValenti fc Fischerl (2005) 
and derived various stellar parameters by matching those 
parameters to theoretical stellar evolu tion models. Their 
values are largely consistent with IValenti fc Fischerl 
(2005): = 0.832 Ig pgg M , = 0.78 ±0.02 R Q , and 




3960 3965 3970 
Wavelength (A) 



3975 



log g = 4.59 



+0.02 
-0.03- 



0.036 - 

We adopt the Takeda et al. values for 



Fig. 1. — Ca II H line for HD 7924. Slight line core emission is 
seen indicating modest chromospheric activity. 

the calculations below. 

Measurements of the cores of the Ca II H & K lines 
show that HD 7924 has modest chromospheric activ- 
ity (Fig. [1]). We measured the chromosp heric activity 
indic es Shk = 0.20 and log i?^ K = -4.89 (|Wright et all 
2004). The full set of Shk measurements for all ob- 
servations of HD 7924 is nearly constant in time (frac- 
tional rms of 1.2%) and does not show a periodicity at 
or near 5.4 d (the planet's orbital period; see §0]). We 
estimate P rot = 38 d an d an age of 3 . 8 Gyr using R^ K 
and B-V ca libration (| Wright et al.ll2004D . Following 
IWrightj (]2005l ). and based on the values of Shk, M Vl 
and B — V, we estimate a radial velocity (RV) jitter of 
2.13 ms -1 . This empirical estimate for jitter accounts 
for RV variability due to rotational modulation of stel- 
lar surface features, stellar pulsation, undetected planets, 
and uncorrected systematic errors in the velocity reduc- 
tion (|Saar et alj Il998t IWrightj [20051 ) . As explained in 
§H1 jitter is added in quadrature to the RV measurement 
uncertainties for Keplerian fitting. 

3. OBSERVATIONS AND DOPPLER REDUCTION 

We obs erved HD 7924 us ing the HIRES echelle spec- 
trometer ()Vogt et al.1 119941 ) on the 10-m Keck I tele- 
scope. The 198 observations span eight years (2001- 
2008) with high-cadence observations — clusters of obser- 
vations on 6-12 consecutive nights — beginning in late 
2005. All observations were made with an iodine cell 
mounted directly in front of the spectrometer entrance 
slit. The dense set of molecular absorption lines im- 
printed on the stellar spectra provide a robust wave- 
length fiducial against which Doppler shifts are mea- 
sured, as well as strong constraints on the shape of the 
spectrome ter instrumental profil e at the time of each ob- 
servation (jMarcv fc ButlerlflOQllValeiiti et al.l[l995Tl . 

We measured the Doppler shift from each star-times- 
iodine spectrum using a mod eling procedure mo dified 
from the method described bv iButler et al. (1996). The 
most significant modification is the way we model the in- 
trinsic stellar spectrum, which serves as a reference point 
for the relative Doppler shift measurements for each ob- 
servat ion. Butler et al. use a version of the iJanssonl 
( 1995) deconvolution algorithm to remove the spectrom- 
eter's instrumental profile from an iodine-free template 
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Fig. 2. — Radial velocity time series for four stable stars in our 
Keck Doppler survey that are similar to HD 7924. These stars 
demonstrate long-term velocity stability and improved Doppler 
performance for observations after the HIRES upgrade in 2004. 
The bottom two panels show only post-upgrade measurements. 
The binned velocities with measurement uncertainties (but not jit- 
ter) are plotted. Panels are labeled with star name, spectral type, 
and reduced rms to a linear fit. 

spectrum. We instead use a new deconvolution algo- 
rithm developed by one of us (J.A.J.) that employs a 
more effective regularization scheme, which results in sig- 
nificantly less noise amplification and improved Doppler 
precision. We defer a more detailed description of this 
method to a forthcoming publication (Johnson et al. 
2009, in prep). 

Figure [5] shows radial velocity time series for four sta- 
ble stars with characteristics similar to HD 7924, demon- 
strating our measurement precision over the past 4-11 
years. In 2004 August, the Keck HIRES spectrom- 
eter was upgraded with a new detector. The pre- 
vious 2K x 2K pixel Tektronix CCD was replaced by 
an array of three 4K x 2K pixel MIT-LL CCDs. The 
new detector produces significantly higher velocity pre- 
cision due to its improved charge transfer efficiency and 
charge diffusion characteristics, smaller pixels (15 /im 
vs. 24 /xm), higher quantum efficiency, increased spec- 
tral coverage, and lower read noise. Our post-upgrade 
measurements exhibit a typical long-term rms scatter 
of ~1.5ms _1 , compared to ~2.5ms _1 for pre-upgrade 
measurements 11 . The pre- and post-upgrade measure- 
ments also lack a common velocity zero point, but we've 
fit for and corrected this offset to within ~2ms _1 for 

11 Measurements prior to JD 2,453,237 are pre-upgrade. 
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Fig. 3. — Time series of binned radial velocities of HD7924 
from Kcck-HIRES with error bars indicating the quadrature sum 
of measurement uncertainties and jitter. 

every star observed at Keck using a large set of stable 
stars with many pre- and post-upgrade observations. To 
further limit the impact of the velocity discontinuity for 
this star, we let the offset float in the Keplerian fits be- 
low, effectively treating pre-upgrade and post-upgrade 
observations as coming from two different telescopes. 

The velocities derived from the 198 observations have a 
median measurement uncertainty of 0.9 ms -1 for single 
measurements. This uncertainty is the weighted stan- 
dard deviation of the mean of the velocity measured from 
each of the -^700 2 A chunks in each echelle spectrum 
([Butler et al.l fl996'). In many cases, we made consecutive 
observations of HD 7924 to reduce the Poisson noise from 
photon statistics. For the Keplerian analysis below (§0]- 
[5]) , the velocities were binned in 2 hr intervals yielding 93 
measurements with an rms of 4.04ms _1 about the mean 
and a median measurement uncertainty of 0.7 m s _1 . The 
binned velocities and associated measurement uncertain- 
ties are listed in Table[4][at the end of this preprint, after 
References] and plotted as a time series in Fig. |3l 

4. ORBITAL SOLUTION 

We searched the radial velocities in Table 0] for the 
best-fit Keplerian orbital solution using the partially- 
linearized, least-squares fitting procedure described in 
IWright fc Howard] (|2008l ). Each velocity measurement 
was assigned a weight constructed from the quadrature 
sum of the measurement uncertainty (listed in Table [4]) 
and a stellar jitter term (2.13ms _1 ). 

The Lomb-Scargle periodogram of the velocities (Fig. 
[4]) shows a strong periodic signal near 5.39 d. We seeded 
the Keplerian search with this period, and a wide variety 
of other trial periods, and found the best-fit orbital solu- 
tion shown in Fig.[5J This minimum in \v — y/xt implies 
a planet in a 5.3978 ± 0.0015 d orbit around HD 7924 with 
a Doppler semi-amplitude of 3.87 ± 0.72 m s _1 and an or- 
bital eccentricity of 0.17 ±0.16. Note that the orbit is 
consistent with circular, consistent with the well-known 
bias against measuring e = in low amplitude systems 
(e.g., lO'Toole et aD l2009f ). The best-fit orbit includes 
a linear trend, dv/dt = —1.07 ± 0.35 ms _1 yr -1 . The ve- 
locity offset between pre- and post-upgrade HIRES ob- 
servations was solved for and is incorporated in Tabled 
Using our adopted stellar mass of 0.832 Mq, we find a 
minimum planet mass of Mp sin? = 9.26 and an or- 
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Fig. 4. — Lomb-Scargle pcriodogram of binned radial velocity 
measurements of HD 7924. The tall peak near P = 5.39 d suggests 
a planet with that orbital period. 



TABLE 2 
Orbital Solution for HD 7924b 



P eir & ii i c t c r 


Value 


P (days) 


5.3978 ± 0.0015 


e 


0.17 ± 0.16 


K (ms -1 ) 


3.87 ± 0.72 


T p (JD - 2,440,000) 


14,727.27 ± 0.87 


u (deg) 


25 ± 60 


dv/dt (ms _1 yr _1 ) 


-1.07±0.35 


Msini (M ) 


9.26 


a (AU) 


0.057 


Nobs (binned) 


93 


Median binned uncertainty (ms -1 ) 


0.7 


Assumed jitter (ms - ) 


2.13 


Reduced rms to fit (ms - - 1 ) 


2.78 




1.21 


FAP (pcriodogram peak) 


<0.0001 


FAP \xu for Keplerian fit) 


<0.001 
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Mass = 9.26 M E /sin i 



P = 5.398 dav 
K = 3.87ms 1 
e = 0.17 




RMS = 2.78 ms- 1 
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Fig. 5. — Single-planet Keplerian model for the radial veloc- 
ities of HD 7924, as measured by Keck-HIRES. The dashed line 
shows the best-fit Keplerian orbital solution representing a 9.26 Mj 
(minimum mass) planet in a 5.398 d orbit. Filled circles repre- 
sent phased and binned radial velocities, while the open circles 
show the same velocities wrapped one orbital phase. The error 
bars show the quadrature sum of measurement uncertainties and 
2.13 m s — 1 jitter. 



bital semi-major axis of a = 0.057 AU. The parameters of 
this orbital solution, and associated statistics, are listed 
in Tabled 

The Keplerian parameter u ncertainties were de rived 
using a Monte Carlo method (|Marcv et al.l l2005bl ) and 
do not account for correlations between parameter er- 
rors. With seven years of observations of a 5 d planet, the 
error on P is quite small (one part in 3600) and largely 
uncorrelated with the errors on other Keplerian param- 
eters. But because K is relatively small, the data are 
nearly compatible with a family of Keplerians ranging 
from circular to slightly eccentric, leading to fractionally 
higher, correlated errors on the other parameters (espe- 
cially e and to). Our estimate for a possible time of mid- 
transit (see §|S]) accounts for these correlated parameters, 
thereby reclaiming some of the clock-like precision of the 
period estimate. 

As we show in §[51 HD 7924 is photometrically stable, 
particularly on the short time scales near the 5.4 d Keple- 
rian period. Thus, the radial velocity signal we detected 
is almost certainly not due to rotational modulation of 
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Fig. 6. — High-cadence radial velocities of HD7924 during a 
long observing run. Measurements were made on 10 of 11 con- 
secutive nights. The velocities show temporal coherence over two 
orbital phases that matches the best-fit Keplerian orbital solution 
to all of the velocities for HD 7924 (the same orbital solution as in 
Fig. [5}. The error bars show the quadrature sum of measurement 
uncertainties and 2.13ms —1 jitter. 

spots or other stellar surface features. 

We considered the possibility that the velocity peri- 
odicity represented by the Keplerian orbital fit in Fig. [5] 
arose from chance fluctuations in the velocities. The am- 
plitude of our orbital solution is less than 4ms~ 1 , which 
is within a factor of two of the quadrature sum of mea- 
surements uncertainties and jitter, so a careful treatment 
is deserved. We tested the null hypothesis in many ways, 
including several estimates of the false alarm probability 
(FAP). 

First, we note that the unreduced value of x 2 decreased 
from 255 for a linear fit to the data (N Q b s ~2 = 91 degrees 
of freedom) to 125 for the single-planet Keplerian fit (85 
degrees of freedom). This corresponds to a decrease in 
from 1.68 to 1.21, and a decrease in the reduced velocity 
rms from 3.81 to 2.78 ms" 1 . 

Figure [5] shows that all of the velocities measured dur- 
ing a recent 11-night observing run appear consistent in 
period, phase, and amplitude with the overplotted best- 
fit Keplerian solution that was derived from the complete 
set of 93 velocities. 

We also computed multiple false alarm probabili- 
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Fig. 7. — FAP determination for HD 7924 using periodogram 
peak height. The histogram shows the peak amplitudes from the 
periodograms of 10 4 synthetic RV data sets, selected with replace- 
ment from the measured RVs. None of the trials yielded a pe- 
riodogram peak amplitude greater than the original fit, implying 
FAP<1(T 4 . 

ties associated with the chance arrangement of ran- 
dom velocities ma s querading as a coherent signal 
(|Marcv et all I2005U iCummingj I2004D . Figure shows 
the periodogram-based FAP test. In this test, we created 
10 4 synthetic data sets by drawing (with replacement) 
the velocities and associated errors from the set of mea- 
sured velocities and using the actual observation times. 
A periodogram was computed for each synthetic data set 
and the peak amplitudes (typically in the range ~5-10) 
are plotted as a histogram in Fig. [7J This distribution is 
cleanly separated from the periodogram peak amplitude 
of the measured velocities, implying a strong coherence 
in the measured velocities well in excess of random fluc- 
tuations. We conclude that it is extremely unlikely (FAP 
<C 10~ 4 ) that a rearrangement of the measured veloci- 
ties would result in a Fourier decomposition with as much 
power in a single mode as we observed. That is, a circu- 
lar or nearly circular orbit of equal or greater amplitude 
would almost never be falsely detected. 

While instructive, the periodogram-based test under- 
estimates the true FAP associated with searching for a 
Keplerian signal because it only measures the strength 
of sinusoidal signals. To overcome this limitation we use 
a test based on \v for a full Keplerian fit, instead of pe- 
riodogram peak height. As with the previous test, the 
synthetic data sets were constructed by drawing with re- 
placement from the measured velocities. Figure [8] shows 
the distribution of \v f° r the synthetic data sets as well 
as the value of \v for the Keplerian fit to the original, 
unscrambled velocities. None of the synthetic data sets 
had a best-fit Keplerian with \v lower than for the un- 
scrambled velocities, implying an FAP for this scenario 
of < 0.001. 

This FAP based on \v is also instructive, but is proba- 
bly too conservative for at least two reasons. (Note that 
these caveats are irrelevant for HD 7924b since we have 
already shown that the FAP is very small, but are im- 
portant for the possible second companion considered in 
§[5] and for future planetary systems.) First, it considers 
all Keplerian solutions equally believable, independent of 
orbital eccentricity. To be sure, nature occasionally pro- 
duces plan etary systems wi th eccentricities of ~0.9 (e.g. 
HD 80606; iNaef et al.ll200l[) that would have seemed un- 
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Fig. 8. — FAP determination for HD 7924 using \ v . The his- 
togram shows Xv from the best-fit Keplerian solution to 1000 syn- 
thetic RV data sets, selected with replacement from the measured 
RVs. None of the trials yielded a \v lower than the original fit, 
implying FAP < 0.001. 

believable not long ago. Yet, when these models arise 
as solutions to synthetic data sets, the high eccentricity 
and low \v are often driven by a few aberrant points that 
sometimes arise in a data set created from non-Gaussian 
noise drawn with replacement. When such solutions are 
proposed for real data, they are often regarded with skep- 
ticism until the eccentric Keplerian hypothesis can be 
tested with additional measurements around the time of 
periastron passage. Thus, in accommodating all Keple- 
rian signals, this test overestimates the FAP. 

The second reason the Xfbased FAP is sometimes 
too conservative is also a failure of over-accommodation. 
This test (as well as the periodogram-based FAP test) as- 
sumes that the observed velocity variation is completely 
due to noise. However we have independent estimates of 
the expected velocity variability of a given star from the 
Doppler analysis pipeline, activity indices, and compar- 
ison stars. For HD 7924, the measurement uncertainties 
(0.7ms -1 ) and jitter (2.13 ms -1 ) predict the level of ve- 
locity rms (2.3 ms -1 ) in the absence of a planetary com- 
panion. Thus, when interpreting this FAP, we should 
discount it by the probability that the observed velocity 
variability is due to noise. 

On balance, we reject the null hypothesis and with 
confidence attribute the observed velocity variation to a 
planet orbiting HD 7924 with a period of 5.398 d. 

5. POSSIBLE SECOND COMPANION 

Motivated by the success of our single-Keplerian fit 
and the higher than expected velocity rms to that 
fit (2.78ms -1 ), we performed a search for a double- 
Keplerian fit. The Lomb-Scargle periodogram of the 
residuals to the single-planet fit (Fig. O shows several 
short-period narrow peaks with modest power (at 12.3, 
16.7, 34.8 d) and an intermediate-period peak with broad 
power (at 147.5 d). Each of these peaks represent a pos- 
sible second planet in the system. 

We seeded the Keplerian search with several trial two- 
planet solutions. The first planet was always seeded with 
the orbital parameters for the best-fit single-planet model 
in §[U The second planet was seeded with the largest 
peaks in the periodogram of single-planet residuals (Fig. 
[§]) and a wide variety of other trial periods. All of the 
parameters were allowed to vary during the fitting pro- 
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Fig. 9. — Lomb-Scargle periodogram of velocity residuals to the 
single-planet fit show in Fig. \E\ The peaks at 12.3, 16.7, 34.8, 
and 147.5 days suggest a possible second planet at one of these 
periods. Additional measurements are needed to tease out such a 
weak signal, if present. 

cess ([Wright fc Howarc3l2008f ). None of these two-planet 
fits changed the orbital parameters for HD 7924b sub- 
stantially. The top four two-planet solutions had best- 
fit second-planet periods of 12.3, 34.8, 16.7, and 147.5 d, 
and best-fit eccentricities of 0.72, 0.76, 0.44, 0.83, respec- 
tively. The high eccentricities raise our suspicion that the 
signals may be due to noise and increase the need for a 
robust tests of the null hypothesis. 

We calculated FAPs for each of these two-planet mod- 
els, building on the methods described in §0J Synthetic 
data sets were constructed by drawing with replacement 
from the residuals to the best-fit single-planet solution 
and adding the coherent, best-fit solution back to the 
scrambled residuals. We performed a thorough search 
for the best-fit two-Keplerian model (as above) on each 
synthetic data set. False alarms were triggered when a 
synthetic data set had a lower value of \ v than the origi- 
nal, unscrambled data. All of the two-planet models had 
an FAP > 20%, even when we restricted the scrambled 
trials to low-eccentricity solutions. Thus, we do not con- 
sider any of these signals to be viable planet candidates 
at the present time. Nevertheless, we will continue to 
hunt for additional planets orbiting HD 7924 by making 
additional observations and by refining our Doppler anal- 
ysis algorithms to reanalyze our extant Doppler spectra 
more precisely. 

6. PHOTOMETRIC OBSERVATIONS 

We acquired photometric observations of HD 7924 with 
the T12 0.8 m automated photometric telescope (APT) 
at Fairborn observatory in southern Arizona. Our bright- 
ness measurements were made between 2006 Decem- 
ber and 2008 October and cover the last part of the 
2006-07, the complete 2007-08, and the first part of 
the 2008-09 observing seasons. The T12 APT and its 
precision ph otometer are very similar to the T8 APT 
described in iHenrvl (|1999t ) . The precision photometer 
uses two temperature-stabilized EMI 9124QB photomul- 
tiplicr tubes to measure photon count rates simultane- 
ously through Stromgren b and y filters. 

The telescope was programmed to gather observations 
of HD 7924 with respect to three nearby comparison 
stars in the following sequence: DARK, A, B, C, D, A, 



SKY A , B, SKY b , C, SKY c , D, SKY d , A, B, C, D. The 
comparison stars A, B, and C are HD4295 (V = 6.39, 
B-V = 0.42, F3 V), HD 10971 (V = 6.94, B-V = 0.49, 
F5), and HD 10259 (V = 6.59, B-V = 1.04, G5), respec- 
tively, while star D is HD 7924 (V = 7.17, B-V = 0.83, 
K0). 

Each complete sequence, referred to as a group ob- 
servation, was reduced to form 3 independent measures 
of each of the 6 differential magnitudes D— A, D— B, 
D-C, C-A, C-B, and B-A. The differential mag- 
nitudes were corrected for differential extinction with 
nightly extinction coefficents and transformed to the 
standard Stromgren system with yearly mean transfor- 
mation coefficients. The three independent measures of 
each differential magnitude were combined, giving one 
mean data point per complete sequence for each of the 6 
differential magnitudes. To filter any observations taken 
under non-photometric conditions, an entire group obser- 
vation was discarded if the standard deviation of any of 
the six mean differential magnitudes exceeded 0.01 mag. 
Finally, we combined the Stromgren b and y differential 
magnitudes into a single (b + y)/2 passband to improve 
the precision. 

Our complete data set consists of 192, 435, and 145 
good group observations from the 2006-07, 2007-08, and 
2008-09 observing seasons, respectively, for a total 772. 
To minimize the effect of any low-level intrinsic vari- 
ation in the three comparison stars, we averaged the 
three D— A, D— B, and D— C differential magnitudes of 
HD 7924 within each group into a single value, represent- 
ing the difference in brightness between HD 7924 and the 
mean of the three comparison stars: D — (A + B + C)/3. 
The standard deviations of these ensemble differential 
magnitudes for each of the three observing seasons are 
1.77, 1.85, and 1.47 mmag, respectively. The standard 
deviation of the complete data set is 1.76 mmag. These 
values are comparable to the typic al precision o f a single 
observation with our 0.8 m APTs (|Henrvill999D . indicat- 
ing there is little or no photometric variability in HD 7924 
or in the three comparison stars. Solar-type stars often 
exhibit brightness variations caused by cool, dark pho- 
tosphcric spots as they a re carried into and out of view 
by stellar rotation (e.g.. iGaidos et ai1l2000f ). However, 
periodogram analysis of the individual three seasons of 
HD 7924 and of the three observing seasons taken to- 
gether, confirms the lack of any photometric periodicity 
between 1 and 100 days. This is consistent with the star's 
modest level of chromospheric activity and its low v sini. 

The 772 ensemble (b + y)/2 differential magnitudes 
of HD7924 are plotted in the top panel of Fig. [101 
Phases are computed from the Keplerian orbital pe- 
riod P = 5.3978 ± 0.0015 d (Table |2J) and the epoch 
JD 2, 454, 727. 99± 0.16, a recent time of mid-transit de- 
rived from the orbital elements. The standard devia- 
tion of a single observation from the mean of the en- 
tire data set is 1.76 mmag. A least-squares sinusoidal 
fit on the orbital period yields a semi-amplitude of only 
0.19 ±0.09 mmag. This very low limit to photometric 
variability on the radial velocity period is strong evi- 
dence that the low-amplitude radial velocity variations 
observed in the star are, in fact, due to planetary reflex 
motion and not to activity-induced intri nsic variations in 
the star itself (e.g. JPaulson et aLll2004f ). 

The photometric observations of HD 7924 near the pre- 
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dieted time of transit are replotted with an expanded 
horizontal scale in the bottom panel of Fig. [101 The 
solid curve shows the predicted time (0.00 phase units) 
and duration (±0.01 phase units) of transits with an ar- 
bitrary drop of 1 mmag from the mean brightness level 
of all the observations. The error bar in the lower right 
of both upper and lower panels represents the mean pre- 
cision of a single observation (1.7 mmag). The horizontal 
error bar immediately below the transit in the top panel 
represents the 2-er uncertainty in the predicted time of 
mid-transit (±0.32 days or ±0.06 phase units). The a 
priori probability of HD 7924b transiting its host star is 
0.06. As discussed in §0 the predicted depth for a transit 
of HD 7924b ranges from 6 mmag down to 0.3 mmag, de- 
pending on the composition of the planet. Transit depths 
of 2-3 mmag or more are easily excluded by the data. 
The 20 observations within the nominal transit window 
between phases 0.99 and 1.01 have a mean brightness 
just 0.06 mmag fainter than the mean light curve level. 
To estimate the limiting depth of detectable transits from 
our current data set, we segregated the observations be- 
twen phase 0.93 and 0.07 into seven bins each 0.002 
phase units in width. We computed the mean depth 
of each bin and compared them to the mean brightness 
level of all the observations. The differences ranged from 
0.46 mmag fainter to 1.12mmag brighter than the overall 
mean brightness. The standard deviation of the seven 
difference values is 0.4 mmag. Thus, we conclude that 
transits deeper than 2-3 mmag do not occur, but addi- 
tional observations around predicted transit times will 
be required to improve the lower limit or confirm transit 
depths around 0.5 mmag or less. 

7. DISCUSSION 

We present the detection of HD 7924b, a super-Earth 
planet with Mp sin i = 9.26 in a P = 5.398 d orbit 
around a K0 dwarf. This paper is the first in a se- 
ries of super-Earth [M sin i < 10 M®) and Neptune-mass 
planet announcements from the NASA-UC Eta-Earth 
Program. The 230 stars in this survey are nearby 
(d<25pc), bright enough for ~lms _1 RV measure- 
ments with Keck-HIRES (V<11; some faint stars suf- 
fer from slightly poorer Doppler precision), have low to 
moderate chromospheric activity (log R' HK < — 4.7), and 
are within two magnitudes of the main sequence. Sev- 
eral candidate super-Earth and Neptune-mass planets 
are emerging from the survey. Some of these require ad- 
ditional follow-up observations to verify the orbit and ac- 
curately estimate Keplerian parameters. Once the survey 
and follow-up observations are complete, the statistics of 
low-mass planets will offer strong constraints on planet- 
formation theory. Just as successful theories of planet 
formation must reproduce the observed planet distribu- 
tion in mass-period space and other p arameterizations 
for J upiter- and Saturn-mass planets (jCumming et all 
2008), they must do the same for low-mass planets. In 
this regard, HD 7924b is an interesting first detection for 
the Eta-Earth survey. Its minimum mass of 9 M® places 
it clearly in the "planet desert", while its small semi- 
major axis (0.057 AU) pushes it into the pileup of short- 
period "hot supe r-Earths" and "h ot Neptunes" seen in 
some simulations (|Ida fc Linll2008l ) . The existence of this 
population depends crucially on an unknown mechanism 
stopping Type I migration and parking the planets in 
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Fig. 10.— Top Panel: The 772 ensemble D - (A + B + C)/3 
photometric observations of HD 7924 in the combined Stromgren 
(6 + 2/)/2 passband, acquired with the T12 0.8m APT over three 
observing seasons and plotted modulo the 5.3978 d orbital period 
of the super-Earth companion. Phase 0.0 corresponds to the pre- 
dicted time of mid transit. A least-squares sinusoidal fit at the 
orbital period yields a semi-amplitude of only 0.19 ± 0.09 mmag. 
Bottom Panel: The photometric observations of HD 7924 near the 
predicted time of transit replotted with an expanded scale on the 
abscissa. The solid curve shows the predicted time of transit with 
a drop in stellar brightness of 1 mmag. The error bar in the lower 
right of both panels represents the mean precision of a single ob- 
servation (1.7 mmag). The error bar immediately below the transit 
in the top panel represents the uncertainty in the predicted time 
of mid-transit (±0.06 phase units). The width of the plot in the 
bottom panel corresponds to this uncertainty in the transit timing. 



short-period orbits. The statistical power of the full Eta- 
Earth survey is needed to check whether HD 7924b (as 
well as other super-Earths listed in Table [2]) is a member 
of this hypothesized class of planets or an outlier in a 
mostly barren patch of planet parameter space. 

HD 7924b is one of only eight super-Earth planets de- 
tected by the radial velocity technique and published (or 
submitted) in the refereed literature (Table[3]). The other 
plane ts are: HD 40307b /c/d, the triple super-Earth sys- 
tem ()Mavor et all [2009); G J 581c/d, the two planets in 
or near the habitable zone (|Udrv et al.l 12001 : GJ876d, 
the super-Earth in the same system as two res onant, 
intermediate-period Jupiters (I Rivera et~aT1 [20051); and 
G J 176b, an 8M ffi super-Earth ([Forveille et al.ll2009f) . 

Table [3] reveals some emerging properties of super- 
Earths, albeit with the limited statistical power of only 
eight detections and through the lens of possible observa- 
tional selection effects. Only three of these planets have 
periods longer than 10 d and only one has a period longer 
than 30 d. All of these planets are in low eccentricity or- 
bits, consistent with circular. As with higher-mass plan- 
ets in short-period orbits, this is probably due to tidal 
circularization. It is curious that HD 7924b and GJ 176b 
are the only Doppler planets in this mass range not in 
known multi-planet systems. This is quite different from 
the 14% of higher-mass pla nets in known multi-planet 
systems (jWright et al.ll2008[ ). 

The host stars of known super-Earths also show emerg- 
ing trends. While Jovian-mass planets are found prefer- 
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TABLE 3 

Planetary and Stellar Parameters for RV-Detected 
Super-Earths (M p isinj < 10 M®) 



Planet 


M p i sin i 

(M(s) 


P 

(d) 


a 

(AU) 


Spec. 




[Fe/H] 


HD 7924b 


9.3 


5.4 


0.057 


KOV 


0.83 


-0.15 


HD 40307b 


4.2 


1.3 


0.047 


K2.5V 


0.77 


-0.31 


HD 40307c 


6.9 


9.6 


0.081 


K2.5V 


0.77 


-0.31 


HD 40307d 


9.2 


20.5 


0.134 


K2.5V 


0.77 


-0.31 


GJ 176b 


8.4 


8.8 


0.066 


M2V 


0.53 


-0.1 


GJ 581c 


5.0 


12.9 


0.073 


M3V 


0.31 


-0.33 


GJ 581d 


7.7 


83.6 


0.25 


M3V 


0.31 


-0.33 


GJ 876d 


5.7 


1.9 


0.021 


M4V 


0.32 


-0.12 



Note. — Planetary parameters (minimum mass, orbital 
period, and semi-major axis) of known super-Earth planets 
(M p i sin i < 10 Mgj) and stellar parameters (spectral type, mass, 
and iron abundance ratio) of their host stars. Only planets 
detected by t he RV metho d are shown. HD 40 307b/c/d were 
discov ered by IMavor et all 1120091) GJ 176 b by IForveille et al.l 
J2009T), GJ581c/d by lUdrv et al.l (I2007T ), and GJ876d by 
IRivera et all p005h . 

entially around higher-mass stars (j Johnson et al.l [2007h . 
the preliminary statistics of super-Earths are consistent 
with the opposite trend. The known super-Earth host 
stars are all K and M dwarfs in the mass range 0.31- 
0.83 Mq. Current Doppler surveys are capable of de- 
tecting super-Earths orbiting G dwarfs (with reduced 
sensitivity), but none have been reported in the lit- 
erature as of 2008 October. Note that super-Earths 
around most F stars and intermediate-mass sub-giants 
with J14>1.3M0 ("retired" A stars) are undetectable 
by the Doppler method due to higher intrinsic RV vari- 
ability observed in these targets. 

The metallicities of super-Earth host stars also show 
the opposite trend seen in higher-mass stars. Note that 
all of the detected super-Earths orbit metal-poor stars, 
in sharp contrast to the positive correlation observed be- 
tween ho st star metallicity and Jupiter-mass planet oc- 
curance ([Fischer fc Valentill2005h . 

Additional planets and non-detections from unbiased 
stellar samples are needed to determine if the above 
trends reflect physical processes or selection effects. 
While planet detectability depends on many factors in- 
cluding the planet's Doppler semi-amplitude K, the am- 
plitude of systematic and astrophysical errors, measure- 
ments precision, the number and timing of observations, 
the number of planets in the system, etc., we can say 
roughly that the lowest-mass detectable planets have K 
comparable to the velocity rms of long-term stable stars 
(see Fig. [2J . Since the Doppler semi-amplitude of super- 
Earths scales as 

3.7 ms- 1 /Py 1/3 /M,\' 2/3 Af pl sin ? 
(1-e 2 ) 1 / 2 \5d) \M & ) 10 M e ' [ ' 

RV surveys are biased against low-mass, long-period 
planets orbiting high-mass stars. For comparision, the 
known super-Earths listed in Table [3] have K = 2.0 
6.5ms- 1 . 

HD 7924 is nearby and bright, making plausible the 
detection of its short-period planet by other techniques. 
While our photometric observations preclude a transit 



down to the level of approximately 0.5mmag during the 
times described in §[51 and the a priori transit prob- 
ability for HD 7924b is 0.06, follow-up campaigns are 
warranted given the extraordinary value of a nearby 
transiting super-Earth. Further, while the radius of 
HD 7924b is unknown, we can place bounds of R p \ — 1.4- 
6 R(f, for planet models s panning pure iron to pure hydro- 
gen ([Seager et al.l 12007). These extremes produce tran- 
sits of depth 0.3 and 6mmag, respectively. A more plau- 
sible silicate- water planet with R p \ ss 2 Rq ([Seager et al.l 
|2007| ) produces a ~0.6mmag transit, near our detection 
limit. 

With an astrometric signature of 0.5 fias, detecting 
HD 7924b astrometrically will only be possible with 
NASA's SIM (Space Interferometry Mission) or its pro- 
posed variants (e.g. Sim-Lite and Planet Hunter) that 
achieve a single measurement precision of < 1 /ias. Ad- 
ditional planet(s) of comparable mass in larger orbits, 
should they exist, would presumably be easier to detect. 

Direct imaging of HD 7924b is infeasible for the fore- 
seeable future given the maximum projected separation 
of ^15 mas for this short-period planet. 

As Doppler velocity surveys probe the observation- 
ally challenging realm of super-Earths and Neptune-mass 
planets, robust detections require continued improve- 
ment in velocity precision, an increased number of ob- 
servations per target, and extra care interpreting the ve- 
locities with careful consideration of the null hypothesis. 
Because the Doppler velocity amplitudes for these lowest 
mass planets are similar in magnitude to the single mea- 
surement precision, we are in an era where it is important 
to demonstrate the fundamental velocity precision and 
long term instrumental stability, and to provide veloc- 
ity data so that others may check the confidence level of 
the claim and explicitly consider the null hypothesis by 
carrying out FAP estimates and other tests. These safe- 
guards will ensure continued confidence and credibility of 
the profound and extraordinary claims of the detection 
of low-mass planets. 
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TABLE 4 
Radial Velocities for HD 7924 





Radial Velocity 


Uncertainty 


JD - 2440000 


(m s~ x ) 


(ms" 1 ) 


12187.970 


5.62 


1.3 


12218.998 


2.86 


1.3 


12220.015 


2.95 


1.7 


12307.770 


9.57 


1.3 


12535.955 


7.59 


1.2 


12575.915 


-2.55 


1.6 


12850.115 


1.19 


1.3 


13239.082 


-2.93 


0.9 


13338.794 


3.20 


0.9 


13400.804 


-1.90 


0.7 


13692.825 


-0.08 


0.5 


13693.006 


-0.09 


0.5 


13693.777 


-0.40 


0.5 


13693.952 


-1.26 


0.5 


13694.820 


0.87 


0.5 


13694.956 


2.40 


0.5 


13695.010 


1.45 


0.6 


13695.898 


5.29 


0.3 


13696.929 


2.99 


0.5 


13724.850 


-4.04 


0.7 


13746.799 


3.51 


0.7 


13748.790 


3.77 


0.6 


13749.785 


5.47 


0.6 


13750.795 


1.05 


0.6 


13751.831 


-1.53 


0.7 


13752.836 


-6.05 


0.7 


13753.801 


-3.46 


0.6 


13775.767 


3.38 


0.7 


13776.805 


7.31 


0.6 


13777.805 


4.78 


0.7 


13778.799 


-0.62 


0.7 
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TABLE 4 — Continued 



Radial Velocity Uncertainty 



J — 2440000 


(ms ) 


(ms 


1 3779 823 


-3.30 


0.7 


13933.124 


6.05 


0.6 


1 3959 1 3D 


0.26 


0.6 


13961.070 


4.64 


0.5 


1 3962 067 


2.93 


0.5 


13963 104 


-3.91 


0.6 


13964.134 


-0.83 


0.9 


1 3982 026 


7.67 


0.5 


1 3983 083 

1 ' ) • .' 1..' ■ ) - V.' O • ) 


5.20 


0.6 


1 3983 987 


-2.50 


0.6 


1 3985 005 


0.89 


0.6 


1 4083 895 


8.07 


0.7 


14084 837 


10.80 


0.7 


141 29 91 7 

Xttx^C^CI I 


-0.30 


0.7 


1 41 30 71 4 


-1.69 


0.7 


141 31 780 


0.76 


0.8 


1 41 38 764 

14:100. 1 Urt 


1.39 


0.7 


14295 1 33 


2.63 


0.7 


14305 135 


3.12 


0.8 


1 4306 1 36 


2.58 


0.7 


14307.122 


-0.45 


0.6 


14308.117 


-6.29 


0.9 


14309 108 


-6.06 


0.6 


14310.108 


-0.86 


0.6 


1431 1 093 

l^yl l.UJU 


1.29 


0.6 


14312.089 


3.27 


0.6 


14313.084 


-3.57 


0.7 


14314 089 


-2.89 


0.7 


14315 106 


1.35 


0.6 


14319.125 


-6.39 


0.5 


14336 084 


-7.84 


0.7 


14337 143 

14UU f • XrtiJ 


-1.27 


1.0 


14339.112 


2.83 


0.9 


14344 047 

X Xv> X ± • VJ'-T I 


2.72 


0.7 


14345 147 

1^'J^Ut Xrt i 


-5.63 


0.8 


14396 864 


-0.96 


1.0 


14397 921 


3.20 


0.9 


14398 932 


1.66 


1.0 


1 43QQ 877 


4.37 


1.0 


14427.831 


2.43 


0.6 


14429 895 


5.23 


0.7 


14430.934 


3.91 


0.8 


14453 738 


-1.64 


0.9 


14454 786 


-4.08 


0.7 


14455 788 


-3.85 


0.6 


14456 771 


0.00 


0.8 


14458.922 


-10.72 


1.3 


14460 768 

x x xu u . i vk.) 


-2.41 


0.9 


14461 695 

x^irtux.ucu 


-0.77 


1.9 


14492.794 


-3.44 


0.5 


14675.111 


-3.67 


0.6 


1 4690 1 1 3 


-5.83 


0.7 


14718.021 


-4.44 


0.6 


14719.063 


-2.22 


0.5 


14720.076 


-0.94 


0.6 


14721.069 


3.44 


0.5 


14722.117 


-1.11 


0.6 


14723.038 


-2.72 


0.6 


14724.006 


-3.75 


0.6 


14725.049 


-0.92 


0.6 


14727.066 


2.59 


0.5 


14727.989 


-0.71 


0.5 



